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Dynamic Kerr effect measurements provide a simple realization of a nonlinear experiment. 
We propose a field-off experiment where an electric field of one or several sinusoidal cycles 
with frequency VL is applied to a sample in thermal equilibrium. Afterwards, the evolution of 
the polarizability is measured. If such an experiment is performed in the Terahertz-range it 
might provide valuable information about the low-frequency dynamics in disordered systems. 
We treat these dynamics in terms of a Brownian oscillator model and calculate the Kerr 
effect response. It is shown that frequency-selective behaviour can be expected. In the inter- 
esting case of underdamped vibrational motion we find that the frequency- dependence of the 
phonon- damping can be determined from the experiment. Also the behaviour of overdamped 
relaxational modes is discussed. For typical glassy materials we estimate the magnitude of 
all relevant quantities, which we believe to be helpful in experimental realizations. 

PACS Numbers: 05.40.jc, 63.50+x, 78.20.Fm 



I. Introduction 

Dynamic Kerr effect studies have been utilized to monitor the reorientational dynamics of 
supercooled liquids for a long time pp. In such experiments, the time-dependent polariz- 
ability is observed after an electric field has been switched on or off ('field on' and 'field off' 
version). In the context of supercooled liquids a comparison of reorientational correlation 
times extracted from Kerr effect data and dielectric relaxation showed that the correlation 
times associated with different Legendre polynomials {1 = 2 for the Kerr effect, / = 1 for 
dielectric relaxation) are very similar 0. 

In a theoretical description of the dynamic Kerr effect, the coupling of the field to 
permanent dipole moments as well as to induced dipole moments has to be taken into 
account. Within a rotational diffusion model, a monoexponential decay ('field off') or a 
biexponential rise ('field on') of the polarizability has been found 0. 

An important property of the Kerr effect is its nonlinearity with respect to the applied 
field. Measurements of nonlinear effects generally provide information about physical prop- 
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erties that is not accessible via linear experiments. An example for this is the distinction 
between dynamic homogeneous and heterogeneous relaxation scenarios A nonlinear 
experiment that was designed particularly in order to make that distinction is nonresonant 
holeburning (NHB) |5j. Here, a large sinusoidal electric field is applied to a sample and the 
linear response of the nonlinearly perturbed system is measured afterwards. The idea is 
to specifically address those dynamical features or relaxational modes which absorb most 
of the energy supplied by the sinusoidal field. This of course depends on the frequency of 
the latter. Whether or not a frequency selective modification of the response is possible 
allows to distinguish between heterogeneous and homogeneous dynamics if the definition 
of reference |4j is used for these terms. 

Realizations of NHB on different systems in most cases lead to the conclusion that 
the dynamics are of a dynamic heterogeneous nature [6j, but also dynamic homogeneous 
behaviour has been reported, e.g. in an ion- conducting glass |7]. All NHB experiments 
performed so far utilized fields with frequencies in the kHz to Hz range, which is ideally 
suited to study slow relaxation phenomena in complex systems like supercooled liquids or 
spin glasses. 

In a recent publication ^ , we presented a model to describe an NHB-experiment on a 
much faster time scale, in the range of inverse Terahertz (THz). These short timescales are 
experimentally accessible since recent advances in the preparation of ultrashort laserpulses 
should allow to provide one or several coherent cycles of a sufficiently strong electric field 
in that frequency range p. We found that in a Brownian oscillator (BO) model frequency- 
selective behaviour can also be expected for vibrational (fast) dynamics. 

Dynamic processes in the THz-range are mainly of a vibrational nature. In glasses, ex- 
cess states compared to the Debye behaviour are found in the vibrational density of states 
(DOS), forming the so-called boson peak. Although there is consensus about the vibra- 
tional character of these dynamics, the lifetime of these states and the damping mechanism 
involved is still a point of interest [Tllj . 

Due to the fact that the Kerr effect signal is of a nonlinear nature, we expect Kerr effect 
measurements (in the THz-range) to be simpler to realize than NHB. Therefore, we study 
the Kerr effect response in this paper, particularly with regard to frequency- selectivity and 
its consequences. In contrast to the purely dissipative diffusion models which describe the 
dynamics on long timescales, we include inertial effects, implement a BO-model |n] for the 
collective THz-dynamics and calculate the Kerr effect response function. The form of the 
signals we find allows to extract the frequency-dependence of the damping, and therefore 
provides important information about the vibrational dynamics. 

The remainder of this paper is organized as follows. In the following section II we 
present the theory. We introduce our model in the first part, followed by a discussion 
about vanishing and nonvanishing contributions to the response in a second part. The 
latter point is deepened in a third part, where we give numerical estimates of the involved 
quantities. Section HI contains the results and a discussion. We summarize and conclude 
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in part IV. 



II. Theory 

1. Nonlinear response functions in the BO-model 

We describe the dynamics of the system in the THz-range by a set of normal modes 
{1} = {<?!; 12, ■ ■ ■ , Qn}- In a molecular system the qi stand for internal degrees of freedom 
(intramolecular modes), but our discussion is focussed on macroscopic systems, especially 
glasses and disordered solids. In this case, the qi are normal coordinates for the collective 
vibrations around the boson peak (intermolecular and intramolecular modes). The classical 
Langevin equation for each mode in the BO-model reads as ^1] 

. , dvi{q}) m,.t ... 

mqi + m-fiqi^ = i^i{t) t: — (1) 

dqi dqi 

Here, 7^ is an Ohmic friction constant associated with the mode i, m is a (mode- 
independent) mass, Ti{t) a Gaussian stochastic force with zero mean and correlation 
i{t)V k{t')) = 26ik'jiP~^m6 {t — f) , = ksT, and V{{q}) the potential energy depending 
on all coordinates. For V we choose the uncoupled anharmonic oscillator potential 

V{{q}) = E^^te) = E^^' {l^" + l^sqf + \e4qt + ...) (2) 

where we scale the anharmonicity constants by the harmonic eigenfrequencies. This 
means that the relative strength of anharmonicities with respect to the square of the corre- 
sponding harmonic eigenfrequencies is constant for all modes. In principle, we could include 
couplings between the different modes in the potential function, but we neglect these effects 
for a qualitative study and consider the system as a set of independent oscillators. 

External electric fields couple to the system via a permanent dipole interaction and a 
polarizability interaction. Due to the fact that the permanent dipole moment is a first 
rank tensor and the polarizability is a second rank tensor, the expressions for the corre- 
sponding energy contributions are proportional to the first and second Legendre polyno- 
mials Pi{cos'd), I = 1,2. Here, is the angle between the direction of the applied field 
and the internal axis of the dipole moments. The energy contributions have the form 
j^perm ^ nEPi{cosi!}) and E^"-'^ = aE'^P2{cosi!}) with the electric field E. The correspond- 
ing external Hamiltonian is thus given by 

Kxt = -/^r [fi{r,{q})E{r,t) P^{cos^) + air, {q})E\r,t) P^icos^)] (3) 

Strictly speaking, n and a in this expression are dipole- and polarizability densities. In the 
following calculations, the applied electric field is chosen as E{y,t) = E * e***^"" sin(f2r) up 
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to tp = 2ttN/Q, and afterwards. Here, denotes the number of cycles applied, E is the 
time-independent field amplitude, and k the corresponding wavevector. The field sequence 
is shown in figure 1. We define the time variable t starting with at tp and we calculate 
the time-dependent evolution of the polarizability after the system has been driven out of 
equilibrium by the sinusoidal field, that is the expectation value 



with the outer brackets (. . .)^ denoting an orientational average. 

We assume that external THz-fields do not affect the orientations i). This is because 
typically the reorientational dynamics in supercooled liquids occur on the timescale from 
microseconds to seconds in the temperature range we are mainly interested in. For com- 
pleteness, however, we also present results for a rotational diffusion model in Appendix 



The (spatially integrated) polarizability in eq.(jl} is calculated in perturbation theory 
with the external field as the perturbation. Its time-dependence is a consequence of the 
interaction with the time- dependent external Hamiltonian (jHl), consisting of two terms de- 
scribing permanent dipole moment and polarizability interactions. As a consequence of 
the fact that orientations are unaffected by the external field the linear response of order 
E vanishes if isotropic systems are considered. This is because the only linear term in E 
is the permanent dipole interaction term in which is proportional to the first Legendre 
polynomial. Together with the second Legendre polynomial in (ji} the contribution van- 
ishes after performing the orientational average (...),? due to the Legendre polynomials' 
orthogonality. 

The lowest nonvanishing order in E is therefore 0{E'^). In contrast to many other 
nonlinear experiments, e.g. NHB, where nonlinear contributions are observed as small 
variations on a usually large linear response background, all the obtained signal in Kerr 
effect measurements is of a nonlinear nature. In detail, the polarizability (jH) in 0{E'^) has 
contributions oc aa from the induced dipole interaction term in (jS)) treated in first order 
perturbation theory, and oc a/i/z from the permanent dipole moment interaction treated in 
second order perturbation theory The orientational dependencies of these terms are 
proportionalities to P2P2 and P2P1P1, both nonvanishing in the orientational average. 

These terms are the contributions found in a rotational diffusion model jS]. In the 
BO-model defined by eq.(0), the interaction with applied electric fields is given by 
Therefore, derivatives of the dipole moments with respect to the replace the dipole 
moments. No external force is acting on the system if permanent and induced dipole 
moments do not depend on the coordinates gj. This need of a coordinate dependency 
of the dipole moments for a coupling of external fields to the system yields a selection 
rule for infrared and Raman activity. For our calculations, we expand dipole moment and 




(4) 
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polarizability in a Taylor series in the qi 

/^(r, {q}) = /^o(r) + E /^'(r)igi + E + ■■■ (5) 

i i 

{g}) = ao(r) + E "'(r)*?* + E ^""(i")^?*^ + • • ■ (6) 

i i 

where primes denote partial derivatives with respect to g^, evaluated at gj = for all 
2 = 1, . . . ,n. For simplicity, we neglect crossterms (oc qiQj) in both expansions. 

The Taylor-expansion of the polarizability (jH)) allows to split the expectation value (jH) 
into a sum Bq + B{t) with 

5o = (P2(cos?9) Jdvaoiv) )^ (7) 

= E^^W = E(^2(cos^)(/dr«:(r) (g,(t)) + i/rfra:'(r) (g.'W) + • • •) ^8) 

Apart from a prefactor, the quantity -Bq + B{t) is the so-called birefringence function. The 
pref actor is given by 27i/{Vn) with the sample volume V and the mean refractive index n 
if a is identified with the difference between the polarizabilities parallel and perpendicular 
to the molecule's symmetry axes [T3] . 

Bq in eq.dTj) does not depend on the normal modes of the system and is therefore 
time-independent in the BO-model. While Bq is the relevant contribution if molecular 
reorientations are considered (as discussed in Appendix B), we have to deal with B{t) 
in the BO-model. According to eq.®, the time-dependent expectation values {qi(t)) and 
(qfit)) have to be be calculated, where the former leads to contributions proportional to a-, 
and the latter to contributions proportional to a'-. Therefore the Fokker-Planck equation 
corresponding to eq.(IH) has to be solved. 

The functional form of the expectation values basically is a consequence of the interac- 
tions with the external field. These interactions are treated in time-dependent perturbation 
theory. In the expression for the external Hamiltonian again the Taylor expansions of 
dipole moment (0) and polarizability ® are used. From these expansions it is clear that 
different derivatives of the dipole moments enter the expressions for the expectation values 
n = 1, 2, . . .. As mentioned above, terms of order E"^ arise from the polarizability 
interaction treated in first order perturbation theory (eq.® below), and from permanent 
dipole moment interaction treated as a second order perturbation feq. fjlOp ). The pertur- 
bation ansatz yields expectation values of the form 

(qnt)) = E'P, E /rfVe'^^'a(^)(r')./f(n,S„t,g (9) 

k=l,2,... •' 



fc,«=l,2,. 
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Here, n = 1,2,..., Pi = Pi{cosi!}) and {a\ , ^\ , /^i ) denote the k-th. or l-th order 
derivatives of a or as in the Taylor expansions ()5|6|1 . Furthermore, fl'{n,Ei,t,tp) and 
fi''^\n, Sj, t, tp) are some time-dependent functions. Here, Sj is an abbreviation for all the 
parameters that characterize the mode, namely uoi, 7^ and the anharmonicities 03,64,.... 
In calculating the functions /f and fl^'''\ the difficulty arises that the solution of the Fokker- 
Planck-equation is only known for harmonic potential. Therefore, the anharmonicities Gj 
must also be treated in (time-independent) perturbation theory. 

The appearance of the different derivatives and the anharmonicities makes the situa- 
tion quite complex. In the following section, we discuss the relevance of these different 
contributions. That means we discuss which derivatives appear in the terms of order aa 
and a/i/i. 

2. Contributions to the Kerr effect response 

In the BO-model, some of the lowest-order contributions vanish. Tables 1 and 2 list the 
lowest nonvanishing contributions for aa- and a/i/i-terms. In this section, we skip the lower 
index i of the derivatives. The terms not listed oc a'a" , Q^a'a' , a'fi'fj,', a'fi"n", a"fi'^", 



contribution 


order 


a'a' 





a"a" 


2 


Bsa'a" 


2 


Qsa'a" 


3 


<d4^a'a' 


2 


e^a'a" 


3 


e^a"a" 


4 



Table 1: Nonvanishing contributions to the response of order aa. A definition of the listed 
'order' is given in the text. 

Q3a"fi'fi', Q^a'fi'fi" and Q^a' ^' ^' all vanish identically in the BO-model, independent of 
the functional form of the applied field. This also holds if one or several of the /i^*^-* are 
replaced by a^'^\ The tables are valid for all Raman experiments or photon echoes, at 
least as long as crossterms in the potential and in expansions of the dipole moments are 
neglected. 

We define a system to count the 'order' of the terms in the tables as follows. The leading 
terms in the forces (derivatives of Tiext with respect to the qi) are n' and a' . These first 
derivatives are defined as order 0, and each additional derivative increases the order by 1. 
Finally, the order of all appearing derivatives is summed up. A similar system is used for 
the anharmonicities. Here, the parabolic potential describes the 'unperturbed' state while 
anharmonicities are considered as perturbations. If a term is proportional to the cubic or 
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contribution 


order 


a fi fi 


i 


a fi fi 


i 


a /i /i 


Q 


Ky^a jj /i 


1 
i 




Q 




Q 


'cy^a n /i 


/I 

4 


Ql.ru" 11' 11' 


Q 


Q^a' jj," ji' 


3 


04aW' 


4 




4 
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Table 2: As in Table 1 for the a/i/x-contributions. 

quartic anharmonicity G3 or G4, we therefore increase the order by 1 or 2, respectively. 
We argue that the lowest order terms in the tables are the most relevant contributions to 
the response, and thus assume that polarizability, dipole moment and potential have fast 
converging Taylor series. In order to justify this approach we have to compare a' to a", /x' 
to yu" and the cubic and quartic anharmonicities G3 and 64 (next section). 

We further assume that terms where different order derivatives of the same quantity 
(a, /i) appear can be neglected. This is because they appear in the form of spatial cor- 
relation functions as e.g. in equation ()10j). Note that this is a common ansatz in glassy 
systems, but the argument is not necessarily valid in molecular systems. Here the deriva- 
tives of dipole moments appear just as numbers, not in the form of correlation functions, 
and there is no general argument for neglecting a term oc a'fi'fi" compared to the term 
oc a"fi'fi', for example. 

3. Estimate of the relevant quantities 
Polarizability and dipole moment 

The Taylor-expansion of the polarizability (jH)) is often truncated after the linear term 
(Placzek-approximation) ^3]. The inclusion of the second derivative frequently appears 
in context with nonlinear Xs-spectroscopy. This is because there is neither a contribution 
oc a' a' a' to the signal (see also 15), nor a contribution oc a'^'fi' (see the table above). 
In an oscillator model, either anharmonicities or the second derivative of the polarizability 
are needed to find a nonvanishing contribution of second order in the permanent dipole 
moment. 

In Kerr effect measurements on liquid water ^H] , observed long lived exponential relax- 
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ations are inexplicable with a purely linear coupling of the field to the g^. From comparisons 
of model calculations with the measured data, the ratio a'^ : a"^ is found to be 50 : 1. 

A simulation of the dynamics in CS2 has been analyzed in an instantaneous normal 
mode approach J7j. The tensorial nature of the polarizability is treated in more detail than 
in our scalar notation. Derivatives of the polarizability are calculated via finite differences. 
For the ratio a'^ : a"^ (in our notation) values ranging from (10 : 1) to (100 : 1) have been 
found. 

We therefore conclude that the second derivative of the polarizability can be neglected 
whenever it is directly compared to the first derivative. Nevertheless, a" has to be taken 
into account for those terms where a replacement of a" with a' yields zero due to the 
properties of the BO-model, see above. 

For the permanent dipole moment, typical values of the first derivative /i' are in the 
range of XDhj A [18j. An estimate of the second derivative /i" is difficult. Here, we simply 
assume that the convergence behaviour of the permanent dipole moment's Taylor series is 
similar as for the polarizability a. 

Anharmonicities 

From a comparison of experimental values on a variety of glasses with the predictions 
of the soft potential model JH] we find a mean value for the quartic anharmonicity of 
G4 ~ 200 A~^. For an estimate of the cubic anharmonicity 63 we use thermal expansion, 
an effect untouched by 64 due to symmetry reasons. In a standard model [201, the average 
thermal displacement of an anharmonic oscillator described by a potential like in eq.(P|) 
is (x(T)) = kBTQ^/ijnujf). We estimate a temperature-dependent lattice constant (or 
atomic distance) a(T) as a(T) ^ a(T = 0) + {x{T)). This leads to 

Aa = a(T2) - a(Ti) = 63^^ (11) 

with AT = T2 — Ti. From experimental data on solid argon [201 (an fcc-crystal) we find 
an expansion Aa ^ O.IA if the temperature is increased from AOK to 80K . The thermal 
expansion is linear in that temperature range. Divided by the average lattice constant, this 
results in a thermal expansion coefficient of ~ 4.6 * 10~"^K^^. A reasonable value for muf 
is O^eV/A"^, also based on JH]- These values with the ansatz in eq. ljllll lead to ©3 ^ QA^^. 
Comparing now the corresponding terms in the equation of motion (0), we find the ratio 

{muj^e^q'^)/{muj^e3q^) = (84 * q)/&3 ~ 30 * q/A. 

Typical values for q are q ~ (0.01 . . .0.1)A, being in the percent range of typical atomic 
distances in solids. Following these arguments, the quartic anharmonicity term is not 
necessarily smaller than the cubic anharmonicity term. Another estimate for O3 based 
on the soft potential model leads to the similar value of B3 ~ 15 A~^. This latter 
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approach is doubtful because the two-level-systems should not be used for an estimate of 
a cubic potential anharmonicity, but all in all we conclude that the quartic anharmonicity 
is not necessarily small compared to the cubic anharmonicity. In our case, we would not 
neglect the term cx Q^a'fi'fi' compared to the term oc Q^a'fi'fi', but the latter term vanishes 
in the BO-model. The next nonvanishing order is oc Q^a" ^' ^' , which we neglect due to 
the relative smallness of a" with respect to a' . It should be mentioned that higher-order 
spectroscopy techniques are more sensitive to higher anharmonicities, as has been pointed 
out in ref.j2I]. 

Summary 

As a result of the proceeding estimates, we are left with the contribution oc a' a' and the 
terms oc a"fi'fi' and oc Q^a'fi'fi'. This finding is analogous to the conclusions from [T3] . 
where the term a'a' is found as the dominating part in the xs signal and corresponding 
terms oc a"a'a' and oc Q^a'a'a' are pointed out as the most relevant contributions to 
the X5-signal. The former term oc a'a' is of minor interest to us, since it represents a 
quasi-linear response, apart from the fact that the field appears quadratically instead of 
linear. In our context we have a squared sine function in the time integral instead of the 
sine function we have in calculating the linear dielectric response (oc /i'yu')- The analytical 
solution for the quasi linear term can be found in Appendix A. The other two contributions, 
denoted as 'harmonic' (oc a"fi'fi') and 'anharmonic' (oc Q^a'fi'fi') terms in the following, 
are discussed in detail in the next section. 

III. Results 

1. Solution for an individual mode 

Inserting the calculated (time-dependent) expectation values {qi) and {qf) into (jHI) we find 
harmonic and anharmonic contributions of the form 

B^'^^"^{t) = E^P,P,P2),j dha"{v),j dh'J dh"f,'{r'),fi'{r"),e^^^'+^"^ /f--(t, tp)(12) 

ganharm^^^ = E^Q^uj^ (P^P^P^) ^ J dha'{v)i j dh' J dh" fi' {r')ifi' {r")ie^^^'''+''"^ /f"'''^"™(t,tp)(13) 

We have set the mass in the Langevin equation to unity. The functions fi{t,tp) are 
linear combinations of (complex) exponentials in time. The exact expressions are given 
in Appendix A. In order to see the general features of the response, it is instructive to 
consider the limit tp — >• cx), when the system has reached a steady state before the field is 
switched off. For the harmonic term, we find 
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where A 1^2 = l/2(7i ± Si = {jf — Aufy/"^. Note that 5i is complex for 7^ < 2uji (un- 
derdamped case, oscillatory motion), and real for 7^ > 2uji (overdamped case, relaxational 
motion). In the case of underdamped motion, the limit tp ^ 00 is achieved if the number 
of apphed cycles fuUfills N ^ while 3> is required in the overdamped case. 



We present the responses in the form of the imaginary part of the time signal's Fourier 
transforms since this representation is best suited to discuss the relevant features. This is 
due to the oscillating form of the responses in the underdamped case. 

If linear response functions are calculated, a simple relation holds between the time- 
dependent response R^^{t) of the system to an oscillating field E{t) = E exp{—iu!t) and 
the Fourier transformed response R^{uj) = dt e*'^*i?^(t) to a pulse field E{t) = E5{t). 
This relation is R^'"(t) = exp{—iu!t)R^{u). It is valid in this form if the AC-field has 
been switched on at time t = — 00. However, since we are dealing with nonlinear response 
functions here, no such relation exists |HI|22]- It is not possible to derive the Kerr effect 
response to the oscillating field from the pulse response. 

Expression (fT3j) shows that the signal's amplitude becomes large if the pump-frequency 
fl is close to the eigenfrequency uJi of the (strongly underdamped) mode. This fact guaran- 
tees frequency-selective behaviour. In addition, signals in the Fourier transform's imaginary 
part (w-dependence) will appear at = 2uJi and 00 = ^^. The singularities in the function 
()15|) stem from the fact that the damping has been neglected partially in deriving this ex- 
pression. The expression for the anharmonic term is more complicated. It is not possible 
to give an expression as simple as eq.()15p for this term, not even in the limit tp — 00. 
Similar as for the harmonic term, signals in the Fourier transforms' imaginary part appear 
at = 2u!i and at = 7^ (due to the terms oc exp(— (Aj + Xj)t) for i = j resp. i 7^ j in 
the time signal, see Appendix A). Furthermore, we find a signal a.t uj = Ui from the terms 
oc exp(— Ajt). Figure 2 shows examples of underdamped mode responses for both harmonic 
(upper panel) and anharmonic (lower panel) contributions. For the harmonic term, the 
discussed limit tp — 00 is shown (dashed line) in addition to results for a single applied 
cycle (full line). For the anharmonic term, we show results of one cycle only for clarity. An 
increasing pump time tp in both cases mainly gives rise to a growth of amplitudes, apart 
from some variations in the relative intensities of the features at u; = Ui, 2ui, 7^. We do 
not discuss in detail the differences between harmonic and anharmonic contributions here, 
although the different functional forms might allow to distinguish between the different 
origins of the nonlinear signal (see also |ISI)- We will focus on the common features in 
the signals. Whether anharmonic or harmonic contributions are of more importance will 
strongly depend on the system examined. 




By further assuming a strongly underdamped mode, terms of order 7j can be neglected in 
the denominators and the Fourier transform's imaginary part becomes 




(15) 
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In addition to the underdamped limit, we consider the overdamped hmit 7^ ^ cjj. This 
hmit can be derived either from eqns. (jl2j) and (jl3p by an expansion in powers of 7"^ or by 
solving the equation of motion of the so-called Ornstein-Uhlenbeck process instead of the 
BO-model. Overdamped motion is relevant also in the Terahertz range, since comparisons 
of experimental data from Raman spectroscopy with theoretical considerations show that 
underdamped as well as overdamped response functions are needed to fit the data, like in 
liquid CHCI3 and CS2 or in CCI4 [2^. We find 

and 

j,anharm,ou^^^^^^ = ^e"(fi) (e"(^^)(l - e-^'^'^f 6-^^'' - 3e"(2^])(l - e-^>*-)e-^>*) (17) 

where e(f2) = Aj/ (Aj — iVt) and Aj = uf /-^i. The appearance of the dielectric loss functions 
again guarantees frequency selective behaviour. Here, the maximum excitation is for 
pump frequency Vt = Ki and Vt = Aj/2 (second term in ^"■"■^"■^"^'OU^^ signals appear at 
frequencies uj = Aj,2Aj. Overdamped mode responses are shown in figure 3 for harmonic 
(upper panel) and anharmonic (lower panel) terms. The curves are for one cycle of the 
pump field. In the overdamped case, the infiuence of additional cycles is very small, 
because the exponentials e"^'*'' in (jl6ll7p are very small after only one cycle if the pump 
frequency is in the range of Aj. The overdamped responses are very similar to results 
from a rotational diffusion model. The resulting expressions for this model are given in 
Appendix B for comparison. 

2. Distribution of modes 

Having solved the problem for the individual modes, we now consider distributions of the 
parameters as a model for macroscopic systems. For the eigenfrequencies Ui a Debye DOS 
g{uji) oc uj'^^m = 2 is a natural choice. We cut off the the distribution at the Debye- 
frequency ojd- Alternatively, a proportionality g{uji) oc uf (m = 4) could be chosen as 
found in the soft-potential model, or a more realistic decay than the abrupt cutoff might 
be used. However, these details have little influence on our results, so we focus on the most 
simple-minded distribution. 

In addition, the damping 7^ in general is frequency-dependent. Among others, inelastic 
X-ray scattering techniques are capable to detect the dependence of sound attenuation on 
momentum transfer Q. Usually, a Q^-law for the damping holds in the region of small 
Q and breaks down at larger values of Q |21] EH] EH] • Together with a linear dispersion 
uji oc Q in the region of small Q we have a quadratic dependence of the damping on the 
eigenfrequencies, '^i{uJi) = buf. Reasonable prefactors are in the range of 0.1 if damping and 
frequency are in meV, as in glycerol b ~ 0.12 meV~^ [5^, in vitreous silica b ^ 0.06 meV~^ 
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or in glassy selenium b ?a 0.17 meV~^ PHI- For glassy selenium, a validity of the relation 
7(Q) oc uJi{QY even beyond the linear dispersion region of Q has been reported PB]. 

Another important point is the Wj-dependence of the appearing dipole moment cor- 
relations in eqs. (jl2ll3|) . Here, we assume the same functional behaviour as found for 
the polarizability correlations, the so-called light-to- vibration-coupling C{uJi). A detailed 
discussion of this ansatz is given in As a consequence we approximately have 

j dha"{v), J dh' J dh"^i\v')i^i'{v")ie'''^^'+^"^ ^ a"^i'^C{u,) (18) 

Different models lead to different results for the function C{ijJi). Typical is a behaviour 
C{ijJi) oc cu" with n = 1 [27j or n = 2 28 . Also noninteger values for n [29^ and a 
nonvanishing limit for cjj ^ (HHI have been reported. 

Using eq. (|T^ . the sum over the different modes i in (jH} becomes the integral 

Bharm^^-^ = j du, a" fx'^ g ico,)C (co,) ft"''"' (t, tp) (19) 

for the harmonic term. We inserted the average value {PiPiP2)-d = 2/15. The anharmonic 
contribution is calculated in the same way by replacing a" with Q-^ojfa' and fi"'^^{t^tp) 
with /f"'^"''™(t, tp) in eq.lfT^. Instead of the time-dependent functions fi(t,tp), the Fourier 
transforms fi{uj,tp) can be used in the integration, leading to the Fourier transformed 
birefringence function B{uj). 

The product g{ijJi)C{uJi) oc tUj-'"^"^ plays the role of an effective DOS in the integral (fT^ . 
Its frequency-dependence governs the weight of different parts of the frequency spectrum. 
The influence of different values is small as long as the sum (m + n) is not too small 
(m + n < 2) or too large (m + ri > 7). In these extreme cases, the frequency dependence 
of the effective DOS is stronger than the functional dependence on a; in the dielectric loss 
appearing in the functions fi{t,tp). That means that only the low- frequency components 
determine the value of the integral for small [m + n), and only high-frequency components 
determine the integral for large {m + n). 

It is important to note that approximated expressions like (jlSp cannot be used in the 
integration (jl9j) . since the neglection of the damping in the denominators produces essential 
singularities in the integrand. 

Figure 4 shows the harmonic (upper panel) and anharmonic (lower panel) contributions 
to the response. We used a Debye DOS, g{uji) oc ojf^uii < ujd = 1 , a light to vibration- 
coupling C{ijJi) = ujf, and a quadratically frequency-dependent damping, 7(0;^) = 0.10;^^. 
For this choice of parameters, all modes are underdamped. One cycle of the pump field 
has been applied, tp = 2-71 /Q (full lines). Dashed curves in the upper panel show results 
in the limit tp — 00, that is an infinite number of pump cycles. The different curves are 
for three different pump frequencies Q. The most apparent difference to the individual 
modes' response is that only one peak is left in contrast to the two (harmonic term) or 
three (anharmonic term) separated signals at = 7j, cjj, 2u;j. This is because the peak at 
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LO = 7j dominates the other contributions after the integration (|19p. The positions u = 2uJi 
(upper panel) and u = Ui (lower panel) are marked by arrows. Only the anharmonic term 
shows small additional features at these frequencies. 

The same observation is made for other choices of the parameters. If the damping is 
linear in eigenfrequencies {'y{uji) = KUi), all modes are underdamped (roughly k < 0.1), 
C{uJi) = uJi and the limit of large pump times is considered (tp oo), an analytical expres- 
sion for the integral (fTIH) can be obtained. Also for this choice of parameters we observe 
that only one peak (at uj = 7^) remains in the spectrum, whereas the other contributions 
become negligibly small. We do not show the result, because the expression is lengthy and 
not instructive, but it shows the correctness of our numerical integrations. The finding that 
the signal appears mainly at the damping is in contrast to our calculations for nonresonant 
holeburning |H], where the main signals appeared at the eigenfrequency. 

We would like to draw the reader's attention to the shifts of the extrema positions 
in figure 4. This is a consequence of frequency-selective behaviour. All modes in the 
chosen distribution can be treated as strongly underdamped modes. As mentioned in the 
individual modes' discussion, most energy is absorbed by oscillators for which the pump 
frequency roughly equals the eigenfrequency, f2 pa cUj. As the signal of these modes appears 
at frequency u = jiiwi) in B^"-^"^{uj) and in b°-^^°''^'^{uj), we find the maximum amplitude 
at ~ 7i(fi), that is ~ 0.1f2^ in the case shown in the figure. These values are marked by 
the dotted vertical lines. The coincidence of the maximum position with 7i(n) holds exactly 
for the curves with infinite pump time tp (dashed lines in the upper panel). The extrema 
in case of a single applied cycle (full lines) are slightly shifted towards smaller frequencies. 
Apart from that small difference, the effect of increasing the number of applied cycles is 
mainly a growth of amplitudes. 

A systematic plot of this observation is shown in figure 5. Here, the observed position 
of the maximum is shown as a function of the pump frequency. Plotted are results for 
a linear as well as a quadratic dependence of the damping as a function of uJi for both, 
harmonic and anharmonic contributions. In agreement with the arguments given above, 
the curves exactly reproduce the assumed frequency- dependence of the damping. 

We finally discuss the separation of the nonlinear a/i/i-contributions from the quasi- 
linear signal oc a'a' . A phase-cycling in order to separate this contribution from the other 
two terms is not possible due to the fact that all terms are of order E"^. This has to be 
contrasted to the situation in NHB where a phase cycle can be devised that allows the 
separation of the nonlinear signal from the total response jn|. For the Kerr effect, we can 
either consider systems having a large dipole moment compared to the polarizability, so 
that the a/iyU-contributions are of more importance than the aa-term. On the other hand, 
in a distribution of eigenmodes, the signal of the aa-term is located at the eigenfrequency 
UJ ~ UJi (underdamped case, see Appendix A) in contrast to the other terms' signal at 
~ 7j. Therefore, the signals should be well separated even if the individual contributions 
cannot be extracted. Note that this latter argument is not valid for overdamped modes. 
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Here, extremum positions of the terms are located very close to each other. 

IV. Conclusions 

We have calculated the Kerr effect response of an ensemble of independent Brownian 
oscillators after a sinusoidal electric field has been apphed to a (glassy) sample in thermal 
equilibrium. Any coupling between the different modes has been neglected, and we assumed 
Ohmic friction for all modes. In principle, it is possible to include coupling effects as well as 
time-dependent damping constants, but we do not expect a qualitative change of our results 
as a consequence of any of these extensions. Our model is predominantly designed for the 
vibrational dynamics around the boson peak in disordered solids, that is the Terahertz- 
range. Permanent dipolc moment interaction (/x) is taken into account as well as induced 
dipole moment interaction (a). 

Prom estimates of the quantities involved we conclude that the dominating terms in 
the signals are proportional to a' a', a"ii'ii' and Q^a'/i'/i' with the cubic anharmonicity ©3. 
Since we are mainly interested in the nonlinear contributions of order a/i/i, we discussed 
the latter two terms, denoted as harmonic and anharmonic terms, in detail. 

All terms were calculated in the classical BO-model, treating the external fields and 
the anharmonicities as perturbations. We considered the Fourier transform of the time- 
dependent responses, more precisely its imaginary part, since this representation turned 
out to be best suited for our discussion. The results contain both hmiting cases, the 
overdamped case and the underdamped case. 

For a single underdamped mode we find signals of similar amplitude at frequencies 
cu — u!i,2u!i,'ji with the modes eigenfrequency cui and its damping constant 7^ for the 
harmonic and anharmonic terms, while peaks occur at a; = Aj,2Aj with Aj = uif/'ji in 
the overdamped case. In addition, we obtain proportionalities to dielectric loss functions 
evaluated at the pump frequency, e"(fi), or twice the pump frequency, €"{2Q), or at least 
functions closely related to the usual dielectric loss. These functions, together with the 
prefactor E'^ (the external field amplitude), substantially represent the energy absorbed by 
the system. It is the appearance of these functions which guarantees frequency-selective 
behaviour if a macroscopic system is considered. 

We model a macroscopic system via a distribution of modes according to a DOS g{uJi). 
We used a Debye model g{uji) oc uof with an abrupt cutoff' at Ui = ujd- The choice of a 
more realistic form for the DOS like the inclusion of a bosonpeak in any form has very 
little influence on our results. Before being able to calculate the response of an ensemble 
of modes, we have to specify additional frequency-dependencies of parameters which are 
irrelevant for the response of an individual mode. These are the frequency-dependence 
of the spatial correlations appearing as prefactors in the response functions of a single 
mode, and the frequency-dependence of the damping, being a parameter in the individual 
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response. The former are determined by the so-called light to vibration coupling C{uJi) in 
case of the polarizability correlations, and we assume the same behaviour for the dipole 
moment correlations. For our model calculations, we used a dependency C{uJi) = uf. As 
above, results for e.g. C{u}i) = Ui hardly differ from the case discussed. For the damping, 
we assumed a proportionality 7i(a;j) oc uj^. Our main finding is that any functional form 
7j(a;j), e.g. r) — 2, can be extracted from experimental Kerr effect data. Due to the 
frequency-selectivity mainly those modes with ~ Q (in the strongly undcrdamped case) 
contribute to the signal. We further find that the signal of a distribution of modes mainly 
shows a peak at the damping {u = '-fi) of these modes, whereas the contributions at 
uj = uji and at u = 2cjj, clearly visible in the individual modes' response, vanish. As a 
consequence, the signal shows a maximum at a; ~ 7i(^^). Therefore, if the time-dependent 
evolution of the system's polarizability after an excitation with a sinusoidal field of pump 
frequency Q is measured for several pump frequencies, the signals are Fourier transformed, 
and the maxima positions are plotted versus Q, the frequency-dependence of the damping, 
7i(a;j), can be determined from the data. We think it would be interesting to compare 
such results to data from other experimental techniques capable to detect the damping's 
frequency dependence, like neutron scattering or inelastic X-ray scattering. 

Although our discussion is focussed on the capability to extract the damping's 
frequency-dependence, our results show that additional information about the dynam- 
ics is contained in the Kerr effect signals, for example information about the origin of 
the nonlinear signal like anharmonicity or higher-order derivatives of the polarizability. 
The principle finding of frequency-selective behaviour allows to specifically address certain 
dynamic features in the sample, which is the main advantage of nonlinear techniques as 
compared to linear spectroscopy. 
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Appendix A 



Full expressions for all contributions 

The exact solutions for the harmonic and anharmonic terms ()12flH|l are 
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(A.2) 



with the eigenvalues Aj defined after eq. fT^ . In the summations, all indices take on values 
1 and 2. 

The term oc a'a' is 



Bf{t) = E'iP.P^), / / dh'a'iT),a'ir%e^^^' ffit,t. 



(A.3) 



with 



ft (^5 



1 



A1-A2 VA2Ai + 4f]2 
which in the overdamped limit becomes 

cql,OU 



1 - e"^2*^)e"^2* 



■e"(2fi)(l-e-^'*'')e 



1 

AlA? + 4fi2 



(A.4) 



(A.5) 



Here, the superscript 'ql' stands for 'quasi- linear'. The numerical value of the orientational 
average is (^2-^2)1? = 1/5. The maximum amplitude of this term is achieved for Vt uji/2 
in the strongly underdamped case and i7 ^ Aj/2 in the overdamped case. 

Fourier transforms fi{uj,tp) for all the functions fi{t,tp) can be derived from the given 
expressions by FT{e'^^) = dte^'^^e"^* = l/(^ — iuj), for ^ = Aj, Ai, A2 or a linear combi- 
nation of the eigenvalues. 
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Appendix B 



Reorientational models 



If reorientational models are considered we can skip the spatial dependence of a in the 
expectation value (^. Here, the number a appears as a prefactor and the influence of 
external fields is a modification of the orientational distribution function (isotropic in the 
equilibrium state). We therefore write the time-dependent polarizability as 



a (P2(cos^9(t)) ) = B''°\t) 



(B.l) 



in contrast to eq.(0]). Note that now the orientations are the time-dependent part. 

The expectation value {P2{'&{t))) is calculated using some dynamic model for the ori- 
entations, e.g. with the rotational diffusion equation 



p{n,t\n') = D 



siwd 



5,9 



+ 



sin 



p{n,t\n') (B.2) 



This equation is the Fokker-Planck-equation for the conditional probability P{VL,t\VL') to 
find the orientation VL at time t, assumed the orientation was VL' at time 0. d and 
denote Eulerian angles, j3 = {kBT)~^ and D is the rotational diffusion constant. The 
external Hamiltonian is again given by eq.(j2I). For the solution of eq. ()B.2|) . the propagator 
P{Q, t\Q') is expanded in a series of Wigner matrices. All details about this expansion and 
helpful orthogonality relations of the Wigner matrices can be found in 

As in the BO-model, the interaction with the external field is calculated using per- 
turbation theory. Note that the corresponding operator in the Fokker-Planck-equation is 
proportional to ^^f^ instead of as in the BO-model. Therefore no derivatives of per- 
manent and induced dipole moment (/x^"-*, a*^"^) appear in the rotational diffusion model. 
Instead, the perturbation terms are proportional to a and /i, and to derivatives of the 
corresponding Legendre polynomial Pi or P2 with respect to i). 

The linear response of order E in B^"^{t) is calculated with the permanent dipole 
interaction of order fiE as a first order perturbation. As for the BO-model, this linear 
response vanishes. The nonvanishing contributions of order E"^ are proportional to aa 
and to ann, where the former corresponds to the quasi-linear response oc a'a' in the BO- 
model, and the latter corresponds to both harmonic and anharmonic terms (see above) in 
the BO-model. We find the Kerr effect reponse 

B^"\t) = l-E'aPe-^^'L^e'\T2,2n)[l-e~^^'^ 



+/?/iV(ri,i^) 3e"(r2-ri,fi) 



(B.3) 
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The dielectric loss is defined as e"{x,y) = ^2^y2 and Vi = l{l + 1)D in the rotational 
diffusion model. Note that the time-dependence is a monoexponential decay with the rate 
r2 = 6D, i.e. the decay rate of the second Legendre polynomial, cf. ref. 0. If we consider 
the limit of many pump cycles or large tp (that means ^ VL/{4:TtD) in this model) we 
find 

B'°Ht, U = oo) = -E2«/5e-^^*e"(3D, ^l) (-a- + -l3iJih"{2D, ^)] (B.4) 
6 \5 D 3 / 

The situation is very similar to the overdamped responses in the BO-model, see 
eqs. (jl6ll7IA.5p . For a description of a macroscopic system, one would choose a distri- 
bution of rotational diffusion constants. Usually, such a distribution is chosen in a manner 
that the linear dielectric response has the stretched exponential form which is often ob- 
served in experiments on glassy materials. As a consequence of the appearance of the 
dielectric loss functions in eq. ()B.3|) ). most of the energy in a distribution is absorbed for 
~ Fi and 2Vl F2, respectively for Vl (2 . . . 3)D. We therefore have for the product 
Titp = l{l + l)D27rN/Q ^ 2'kN for the mainly energy absorbing modes. Therefore, the 
appearing exponentials exp(— F^tp) in eq. ()B.3|l are practically zero after only one applied 
cycle, meaning that the limit ()B.4|) is always a good approximation, very similar as for the 
overdamped BO-model response. 

Instead of rotational diffusion, other reorientational models can also be considered, for 
example isotropic or anisotropic random jump models. We expect a similar form of the 
Kerr effect response also for these models. 
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Figure captions 

Fig.l : Electric field sequence: One or several full cycles of the form E{t) — Esin (Qr) 
excite a sample in thermal equilibrium. In the text, Q is denoted as 'pump frequency' 
and tp = 2NTT/Q as 'pump time' with denoting the number of cycles. After the field 
is switched off at tp, the time-dependent evolution of the polarizability is monitored. 

Fig. 2 : Imaginary part of the Fourier transformed response of a single underdamped 
mode for the harmonic term (upper panel) and anharmonic term (lower panel) . The 
eigenfrequency is u;^ = 1 and the damping 7j = 0.1. The pump frequency is Q = 1. 
For the harmonic term responses of one (full line) and an infinite number (dashed 
line) of applied cycles are shown. Because the infiuence of increasing the pump time 
is mainly a growth of amplitudes we show only the response of one applied cycle 
for the anharmonic term (lower panel). The features' positions cu — ^i,uji,2uji are 
marked by vertical lines. 

Fig.3 : As in figure 2 for an overdamped mode. Here is a;i = 1 and 7^ = 10, leading to 
a Ornstein-Uhlenbeck rate of Aj = 0.1. The pump frequency is Q = 0.1. One cycle 
has been applied. 

Fig.4 : Upper panel: Imaginary part of the Fourier transformed signal [i?''™™(tc;)]" from 
a distribution of modes (harmonic term). The mode distribution function is a (nor- 
malized) Debye, g{uJi) oc cufjCUi < 1. The light to vibration coupling is C{u!i) — cuf, 
and the damping is chosen as ^i{uji) — 0.1 uf. Shown are results for three different 
pump frequencies fl = 0.01,0.03,0.1 from left to right, in each case for one applied 
cycle, A^ = 1, (full lines), and for an infinite number of cycles, N ^ 00 (dashed lines). 
The frequencies uj = 7i(0) are marked by vertical dotted lines (extrema positions of 
the corresponding responses), and marked by arrows at a; = 2fl. No features are 
visible at the latter positions. The pref actors 2/15 E'^o;"//'^ were set to unity. 
Lower panel: As in the upper panel for the anharmonic term, [S""'*"^"*(a;)]". Here, 
minor features are observable at cu — (marked by arrows) and at cu — 2 ft. Only re- 
sults for one cycle are shown. Note that amplitudes are plotted negative. Prefactors 
set to unity are 2/15 E'^Q^a'/i''^. 

Fig.5 : The extrema positions of [B'^"'^"^{u)]" and [S""'^"''"*(a;)]", plotted versus pump 
frequency Q. DOS and fight to vibration coupfing are chosen as in figure 4. The 
upper line is for frequency-dependent damping '~fi{uji) — 0.1 Ui, the lower lines for 
7j(c<jj) = 0.1 ujf. The functions exactly reproduce the assumed frequency-dependence 
of the damping. Results of harmonic and anharmonic terms are not distinguishable. 
Full lines are for one cycle of the pump field. The dashed line stems from the harmonic 
term in the limit of large pump times, illustrating the slight shift towards larger uj 
with increasing tp mentioned in the discussion of figure 4. 
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